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Chiral limit of QCD 

Rajan Gupta ^ 

''T-S Group, MS B285, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 U. S. A. 

This talk contains an analysis of quenched chiral perturbation theory and its consequences. The chiral behavior 
of a number of quantities such as the pion mass mj, the Bernard-Golterman ratios i? and x- tlif^ masses of nucleons, 
and the kaon B-parameter are examined to see if the singular terms induced by the additional Goldstone boson,?;', 
are visible in present data. The overall conclusion (different from what I presented at the lattice meeting) of this 
analysis is that, with some caveats on the extra terms induced by ?/ loops, the standard expressions break down 
when extrapolating the quenched data with niq < mg/2 to physical light quarks. I then show that due to the 
single and double poles in the quenched ?/, the axial charge of the proton cannot be calculated using the Adler- 
Bell-Jackiw anomaly condition. I conclude with a review of the status of the calculation of light quark masses 
from lattice QCD. 



1. INTRODUCTION 



The main question this review attempts to 
answer is "should the ostrich care about the 
alarmists view of quenched QCD"? The alarmists 
are two groups, Sliarpc. Labrenz, and Zhang [3] 
[16] [18] and Bernard and Golterman [1] [2]. They 
have calculated, using quenched chiral perturba- 
tion theory, a number of quantities to 1-loop and 
find that in the quenched approximation rf loops 
give rise to unphysical terms in the chiral expan- 
sion and that in many cases the chiral limit is sin- 
gular. Also, the coefficients in the chiral expan- 
sion (including those of the normal chiral logs) 
are different in the full and quenched theories. 
The ostrich is everyone who wishes to continue 
using the chiral expansions derived for the real 
world for extrapolating quenched data to the chi- 
ral limit. The answer is, unfortunately, YES they 
should care. 

The artifacts due to rj' loops can potentially 
invalidate all the extrapolations to the chiral 
limit. The hope is that since these are loop cor- 
rections and potentially large only in the limit 
TJiq — > 0, therefore, there might exist a window 
in niq where the leading order chiral expansion 
is valid and sufficient, albeit with coefficients dif- 
ferent from those in full QCD. Extrapolations of 
the quenched data from this range to the phys- 
ical light rriu may prove to be sensible, and the 



difference between the full and quenched coeffi- 
cients taken as a measure of the goodness of the 
quenched approximation. With this goal in mind 
I analyze the existing quenched data in the range 
m,,/4 — /Ji,, and show that terms induced by the 
rj' are already visible and statistically significant. 

In Section 9 I switch gears and review the sta- 
tus of m and /Ji.,. The quenched Wilson fermion 
data for m is almost a factor of two larger, even 
at /3 = 6.4, than that for quenched staggered or 
nf = 2 staggered or Wilson fermion data. The 
estimates of rrig depend on whether K or K" or (j) 
is used to set the strange scale. These systematic 
differences are much larger than statistical errors 
and need to be brought under control. 

2. QUENCHED CHIRAL PERTURBA- 
TION THEORY 

Morel [5] gave a Lagrangian description of 
the quenched theory by introducing ghost quark 
fields with Bose statistics. This Lagrangian ap- 
proach has been further developed by Bernard- 
Golterman into a calculational scheme. To the 
order we will be concerned with jCbg is 



— str 



+ aodf,<^odf,<^o - ml<^l (1) 

where / = = 131 MeV is the pion decay 
constant, S = exp(2«n//), M is the hermitian 
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Figure 1. The pseudoscalar propagator, (b) the 
hairpin vertex, and (c) the one bubble contribu- 
tion to tlio ■// propagator in full QCD which after 
summation of all diagrams has the form shown. 
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quark mass matrix, fi sets the scale of the mass 
term, and str is the supertrace over quarks and 
ghost quarks. The last two terms involve the field 
$0 = — fi')/'\/2, where ?)' is the ghost (com- 
muting spin- 1 / 2 ) field companion to the 77' . These 
terms are treated as interactions and give rise to 
"hairpin" vertices (see Fig. 1) in the rj' propa- 
gator. This introduces two new parameters, m,^ 
and a momentum dependent coupling oojp. in 
the quenched analysis. In the full theory this ver- 
tex and the tower generated by the insertion of 
bubble diagrams sum to give 77' its large mass, 
711^/(1 — ao), while in the quenched theory the 
rj' remains a Goldstone boson and its propagator 
has a single and double pole. 

The strength of the vertex, -mg, has been cal- 
culated on the lattice by the Tsukuba Collabora- 
tion [4] by taking the ratio of the disconnected to 
connected diagrams. It has also been determined 
using its relation to the topological susceptibility 



^0 = 2nfXt/f^ 
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measured on pure gauge configurations. These 
methods give 750 < mg < 1150 MeV. The pa- 
rameter that occurs repeatedly in the chiral ex- 
pansion of quenched quantities is (5 = ■/7ig/247r^/J. 
Using = 131 MeV and the above estimates for 
mo gives 0.14 < (S < 0.33. Current quenched data 
supports a value between 0.1 < S < 0.15; differ- 
ent lattice observables give varying estimates due 
to statistical and systematic errors. 

Let me first give an intuitive picture of why 
the rj' propagator gives extra contributions. The 
enhanced logs due to the rj' are infrared diver- 
gent, so it suffices to consider the = limit in 



the rj' propagator. The single pole term is akin 
to the pion in the full theory, 1/m^. while the 
double pole term (due to the hairpin vertex di- 
agram) is -^ma^^. Thus any time there is a 
normal correction term like m^Lnm^ from pion 
loops there will also be a singular term of the 



form 



;Lnr 



m-gLnm.^ ^ r^Lnm^. This is 



exactly what one finds in the chiral expansion for 
m^. Similarly, in the case of rrinucleon the reg- 
ular chiral correction is c< m,^, and the rf gives 
an extra term on mgm^. My goal is to expose 
these extra terms in the present lattice data for 
different observables, and extract ^ from them. 

Further details on the formulation of the 
quenched chiral lagrangian and on the cal- 
culation of 1-loop corrections are given in 
Refs. [1] [3] [16]. The 1-loop corrections in 
the full and quenched theories show that 

• the expansion coefficients are different, 

• there are enhanced chiral logs, 

• there are no kaon loops with strange sea quarks, 

• values for parameters like /. //,, .. are different 
in the quenched expressions. I will assume that 
this difference is implicit in all subsequent dis- 
cussion even when the same symbols are used for 
the two theories. Before addressing the conse- 
quences of these differences for the various physi- 
cal quantities and their significance in the present 
data, I would like to mention the difference in the 
strategies, after 1-loop corrections have been cal- 
culated, of the two groups of alarmists. I find 
that knowing their respective emphasis helps in 
reading their papers. 

Sharpe and collaborators focus on determining 
quantities that can be extracted reliably from 
quenched simulations. Using real world values to 
determine the chiral parameters (or commonly 
accepted ones if these are unknown parameters in 
XPT) they require that the chiral corrections are 
small in both the full and quenched expressions, 
as well as in their difference. Observables satisfy- 
ing these conditions are the "good" candidates. 
Bernard and Golterman concentrate on testing 
quenched xPT by forming ratios of quantities 
which are (a) free of O(p^) terms in Ccpt and 
(b) independent of the ultraviolet cutoff used to 
regularize loop integration. The quenched chiral 
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expansion of such ratios have terms proportional 
to the extra parameter S. Since these terms can 
be singular in the chiral limit, it is necessary to 
assume that there exists a window in rrig where 
the 1-loop result is reliable. Then 5 can be de- 
termined from fits to the quenched expression 
provided the fits to the quenched and full theory 
are significantly different. 

3. mi VERSUS m„ 

TV 'I 

Gasser and Leutwyler [6] [11] show that in full 
QCD 

m,l = (1 + ^L{m^) - + 0(m,)) (3) 

where L(m) = m^Ln(m^/A^)/87r^/^. Bernard 
and Golterman [1] and Sharpc [3] show that these 
logs are absent in the quenched approximation. 
Instead, for ao = 0, they get 

2 

K)q = 2^m,{l-Shn(^) + ...). (4) 

whore A is some typical scale of xSB. This ex- 
pression has been refined by Sharpe, who summed 
up the leading logs for the degenerate case m„ = 
rrid = rrig. We use his result [10] 

j^j^ (™7r)Q = Co - —^—Lnm„ + cim„ + C2m^ (5) 
nig 1 + d "J ' 

to extract 8 from a compendium of staggered 
fermion data at /3 = 6.0 [7] [8] [9]. Expressing all 
quantities in lattice units, the best fit gives 

Ln^^^^ = 1.54-0.044Lnm<,-hl.2m<,-2.8m^(6) 

This implies that 8 « 0.053, i.e. much smaller 
than the value « 0.3 based on full QCD param- 
eters; however, the breaking of fiavor symmetry 
in staggered fermions has an interesting conse- 
quence for this analysis. The // operator is a 
singlet under staggered flavor, and different from 
the Goldstone pion which has fiavor 75 . Thus one 
should use the corresponding non-Goldstone pion 
mass in terms that come from the //. In Fig. 2 the 
fit uses the tt (which has fiavor 7475) mass in the 
log term as it is better measured and consistent 
with the fiavor singlet case. The result is 

Lj^II^Oq ^ 1.35-0.13Lnm|-hl.5m|-2.4m1^.(7) 



Figure 2. Fit to staggered m,^/mg data. At /3 = 
6.0, an estimate of strange quark mass is m^a « 
0.025 
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In this form the coefficient of the Lnm^ term 
is 8. Thus 8 « 0.13, a value consistent with 
the estimate 0.14 based on the calculation -mo = 
750 MeV. Also note that since the mass of the 
flavor singlet state, tt, does not vanish as rriq 0, 
therefore, there is no singularity at finite a due to 
the enhanced logs. 

The above analysis shows that if one wanted to 
extract the value of in the expansion = 
j4^mq -|- . . ., then the quenched data would give 
a significantly different result depending on the 
kind of fit used. If one assumes that the 5 data 
points by the Staggered collaboration [7] repre- 
sent a window in which is valid and chiral 
corrections are negligible, i.e. the relation = 
j4^mq is sufficient (as expected at small enough 
rUg in full QCD), then one gets r?i^ = 5.87mq [7], 
whereas Eq.7 gives ~ 3.9, a significantly dif- 
ferent value. The fit in Eq. 7 shows that over a 
range of m^, the chiral log and higher order terms 
can conspire to produce a fiat region. 

Finite size effects in increase the value of 
(m^)Q/mq. so one might attribute the 4% devi- 
ation at nig = 0.0025 in Fig. 2 to this artifact. 
Fortunately, Kim and Sinclair [8] have obtained 
high statistics data for nig = 0.0025, 0.005, 0.01 
on lattices of size L = 16, 24, and 32 as shown 
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in Fig. 2. There is clear indication of finite size 
effects on L = f6 lattices, but the near agree- 
ment between X = 24 and 32 data confirms that 
L = 32 is essentially infinite volume. To conclude, 
the data show that the lowest order chiral expan- 
sion has broken down and the effects of // logs 
are manifest for rrig < mgj2. A similar analysis 
with Wilson fermions is not yet useful because 
the lowest niq used in simulations is <^ 0.4m,,, i.e. 
the point where staggered fermions just start to 
show significant deviations. 

4. BERNARD- GOLTERM AN RATIO R 
AND U 

The chiral behavior of in full QCD has been 
analyzed by Gasser and Leutwyler [11] to be 



u = 



1 

-L{mK) + 



/(m„ + + TO5)I/4 + m^I/s] (8) 

where I/4 and I/5 are two 0{p^) constants they 
introduce. In the quenched theory, with ctQ = 0, 
Bernard-Golterman and Sharpe get 

U = /(l + m^Ls). (9) 

The absence of pion and kaon chiral logs in the 
quenched expression is a 13 — 19% effect (corre- 
sponding to the range A = 0.77 — 1 GeV for the 
chiral symmetry breaking scale in L(m)) using 
full QCD parameters! To reliably compare full 
and quenched theories without ambiguities due to 
the cutoff A and O(p^) terms. Bernard-Golterman 
construct, in a 4-fiavor theory, the ratio 



R 



(10) 



/ll'/22' 

where mi = mi; and 1112 = m2<. The xPT ex- 
pression for i? in the full and quenched theories is 



R^ = l + 
R"^ 



1 



32772/2 
l + S 



,Ln- 



-Ln- 



+ m22'Ln- 



(11) 



'"-11' "'■22') ""22' 

where the quantity within [ ] (called X) increases 
with the mass difference m2 — mi. 

The quenched Wilson data for R obtained by 
the LANL [12], UKQCD[13], and Bernard et 



Figure 3. The Bernard-Golterman ratio R versus 
the full QCD expression given in Eq.ll. 
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a/. [14] collaborations are shown in Figs. 3 and 4 
versus the full and quenched expressions given in 
Eq.ll. The slope of the fit to i2<3 gives 8, while for 
R^ the expected slope is unity. The data favor 
the quenched expression and give 6 = 0.10(3). 

The caveat in this case is that the two points 
at largest X*^ are obtained with m,2 = I'm^. so 
one could argue that 1-loop xPT is not reliable 
for these masses. Barring this technicality, I 
believe that this quantity provides the cleanest 
determination of 8. 

5. BERNARD-GOLTERMAN RATIO x 

AND (V"/') 

Bernard-Golterman construct a second quan- 
tity that is independent of A and O(p^) terms 

{dd) (Ml.-Ml^\ {-ss) ^^2^ 



(13) 



^ {uu) V J {uu) 

for which y^T gives 
rris - md 



Xtr 



r niu nid - m„ m^ 

XQ = Xtree + 8 Lu Lu 

' nid rris - m„ m^ 
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Figure 4. The Bernard-Golterman ratio R versus 
the quenched expression given in Eq.ll. 
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To evaluate these expressions requires data for 
the condensate at three values of rriq and pseu- 
doscalar masses for the combinations tt = uu, 
= sd. = .S7t. At present only the stag- 
gered [7] and Wilson [12] fermion simulations at 
P = 6.0 by the LANL collaboration have all the 
necessary data. Our results for (5 = (x — )(tree)/Y, 
where Y is the factor multiplying S in the expres- 
sion for xq in Eq.l4. are given in Table 1. 

The staggered data have large errors and would 
give the wrong sign for 5. (I have not taken 
into account the difference between Goldstone 
and non-Goldstone mass in terms that come from 
rj' loops.) With Wilson fermions the condensate 
in the chiral limit can be calculated in two ways, 
using the GMOR relation or the Ward Identity 
as explained in Ref. [15]. At finite rUg there are 
lattice artifacts which we cannot control, never- 
theless, the data give reasonable value for S. This 
is probably fortuitous and I believe that much 
better data is needed in order to extract S from 
the chiral condensate. 



Table 1 

The Bernard-Golterman ratio X 

Staggered Wilson(GMOR,) Wilson(WI) 



X 


0.549(30) 


0.608(6) 


0.614(5) 


^tree 


0.517(14) 


0.620(2) 


0.620(2) 


XF 


0.509(15) 


0.616(2) 


0.616(2) 


s 




0.10(5) 


0.05(4) 



6. CHIRAL EXTRAPOLATION OF THE 
NUCLEON MASS 

The behavior of baryon masses has been calcu- 
lated in xPT and has the general form [6] 



+ 0(m^Ln'm.n-) 



(14) 



where Mi are tt, K, rj meson masses. The term 
proportional to M^ comes from pion loops and 
is 25% - 50% of Mb for the octet. For exam- 
ple, using the results of Bernard et al. [19] one 
finds Mjv = 0.97 + 0.24 - 0.27 respectively for the 
first three terms in Eq.l4. Thus, the loop cor- 
rections in individual masses are large and one 
could question whether xPT is applicable at all 
to baryons. On the other hand xPT results for 
mass differences and the Gellmann-Okubo for- 
mula work very well, just as in the quark model. 
So, it is possible that the loop effects somehow 
conspire to just shift the overall scale, in which 
case xPT is useful and it is worthwhile examining 
the consequences of the quenched approximation. 

Labrenz and Sharpe [16] have extended the 
Lagrangian approach of Bernard-Golterman to 
baryons using the "heavy-quark" formalism of 
Jenkins and Manohar [17]. They show that along 
with a modification of the Ci in Eq.l4 one gets 
a mgm^ term due to rj' loops. The quenched ex- 
pression for degenerate quark masses is (assuming 
Oq = 7 = 0. where 7 is a parameter in the baryon 
sector of Ccpt and defined in [16]) 

Mb=M + c^^hM^ + c^'^^Ml + c^^^Ml + ... (15) 

where cM'^ 2.5. r^^) ~ 3.4. and c^^) ~ -1.5 

using full QCD values for the parameters. 

Fits to lattice data using Eq. 15 are not very 
reliable because the number of light quark masses 
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Figure 5. Fit to the LANL nucleoli mass data. 
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explored are typically 3 — 4 and only the point at 
the heaviest mass (typically > 2m,,) shows 
any significant deviation from linearity. I find 
such 4-parameter fits to 4 points very unstable. 
For example, in the case of LANL data [12]. oven 
the different JK samples give completely different 
values of c^*' . The best I could do was to fix one of 
the parameters and make a 3-parameter fit and 
then vary the fixed parameter to minimize x^- 
The best fit (obtained by fixing any one of the 
less well determined coefficients, cMK c^^' or rS^\ 
as one gets the same final result on minimizing 
X^) to the LANL data expressed in units of GeV 
is shown in Fig. 5 and gives 

Mb = 1.16-0.36M^ + 1.6M2-0.5M^ (16) 

Assuming fM"^ = —2.5. Eq.l6 gives S « 0.14. The 
same method applied to "012 sink'" data from the 
GFll collaboration[20] at /3 = 5.93 gives ( up- 
dated version of the fit presented in Ref.[16]) 



Mb = 1.18 - l.OM^ + 3.0M^ - 1.3M^ 



(17) 



which implies that S « 0.4, and c^^' and c^^' have 
values close to those for full QCD. 

7. THE KAON B PARAMETER 

The kaon B parameter is a measure of the 
strong interaction corrections to the if" — 



mixing. It is one of the best measured lattice 
quantities. For details of the phenomenology 
and of the lattice methodology I refer you to 
R,efs.[21][22][23]. Here, I present a summary of 
just the chiral behavior. 

Zhang and Sharpe [18] have calculated the chi- 
ral behavior of Bjc in both the full and quenched 
theories . The full QCD result is [23] 



Bk = B 



l-i3+j)yLny+by+ce^y+0{y^] 



(18) 



where y = mj^/^Sir f ) « 0.2 and e = (m^ — 
mj)/(m,^ + m,i) measures the degeneracy of .s and 
d quarks. B is the leading order value for Bk, 
which is an input parameter in xPT, and b and c 
are unknown constants. The quenched result [18] 



B 



K 



bQ 



l-(3+e2)yLnj/+6S + c' 



+5 



2e 



1 + e 



+ 2 



(19) 



has exactly the same form except for the addi- 
tional term proportional to 5, which is an artifact 
of quenching. The term proportional to 8 is singu- 
lar in the limit e — ^ 1, therefore extrapolations of 
quenched results to the physical non-degenerate 
case are not reliable. For e = this term vanishes, 
so unless one works close to e — ^ 1 (for which there 
is little incentive in the quenched approximation), 
it is unlikely that we will, in the foreseeable fu- 
ture, be able to extract 5 using Eq.l9. 

The constants B, b, c are different in the full 
and quenched theories and cannot be fixed by 
XPT. Assuming B = B'^, the coefficient of the 
chiral log term is the same for f = 0. This is the 
best agreement one can expect between the two 
theories. As a result Sharpe [23] advocates that 
Bk with degenerate quarks is possibly a ''good"' 
quantity to calculate using the quenched theory, 
though systematic errors due to use of degenerate 
quarks are hard to estimate. 

Using full QCD values, 3yLny « 1, so one can 
ask whether this normal chiral log is visible in the 
present data and whether it should be included 
in the extraction of Bk^ With existing data it is 
hard to distinguish this term from the one linear 
in y as the range of mfc is not large enough to 
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Figure 6. Evidence of finite size effects in en- 
hanced chiral logs in By- 
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significantly affect the logarithm. Also, there ex- 
ist data for mq « ■m.,,/2, so for degenerate quarks 
(which, as explained above, is the best one can 
do with the quenched theory) there is no need for 
an extrapolation. 

For staggered fermions B^: can be written as 
the sum of two terms, Bk = By + Ba-, each of 
which can be analyzed using xPT. These quanti- 
ties are defined in Ref. [3] and are explicitly con- 
structed such that they do not diverge as Ijm^- 
in the chiral limit. Both By and Ba have en- 
hanced logs (terms proportional to Ln?/ and not 
suppressed by powers of y) that have nothing to 
do with quenching, i.e. are not due to the rj' . It 
is these logs, or more precisely the volume depen- 
dence of these logs, that has been seen in lattice 
data. Sharpe [3] has shown that this volume de- 
pendence is of the form 

Bv{L)-Bv{^) = -{Ba{L)-Ba{^)) (20) 

rriKL Stv'^P ' 

The constant &2 is not well determined, but the 
shape of the rrifc dependence is. The staggered 
fermion data at /3 = 6.0 on 16^ and 24^ lattices 
[24] are shown in Fig. 6 and qualitatively confirm 
the expected finite size effects in the chiral logs. 



8. MATRIX ELEMENT OF SINGLET 
AXIAL CURRENT IN THE PROTON 

Ever since the measurement of the spin struc- 
ture of protons using deep inelastic muon scatter- 
ing from protons by the EMC collaboration[25]. 
there has been much interest in the calculation of 
the forward matrix elements of the singlet axial 
current in the proton, (p. .s|(7v'7^,75q|p. .s). There 
are two possible Wick contractions that con- 
tribute to this matrix element (ME). These con- 
nected and disconnected diagrams are discussed 
in [27]. Since the disconnected diagram is hard 
to measure, Mandula [26] used the anomaly con- 
dition to derive the relation 



(J},s\Ai^\p,s)si^ 



Nf^ lim 



2tt if-s-o q ■ s 
(j>',s\TrF^,F^,iq}\p,s) (21) 

where q = p—p' and s is the proton's spin vector. 
The hope was that it would be easier to measure 
the ME of this purely gluonic operator. Since 
the rj' propagator contributes to this ME at tree 
level, the question arises whether Eq.2f is valid 
in the quenched approximation. The answer is 
NO [27]. Consider the Fourier transform of the 
anomaly relation 

kfi {?', s\ Af,{q) \p. s) = 2m, (p". s\ P \p. s) + 

Nf^ {p',s\TrFF\p,s). (22) 

Each of the three ME in Eq.22 can be parame- 
terized in terms of form-factors as 

{p'. .s| Aij{q) \p,s) = m7^75uG^ - iq^u^y^uGf, 

{p',s\P\p,s) =U^5UG^, 
{p', s\Tr FF \p, s) =u-f5uG^. (23) 

In the quenched approximation the singularities 
in these form factors for on-sliell ME with respect 
to q"^ and due to the // propagators are 

G^{q'^) no rj' poles, 

^A/ 2\ «2 , fli 



Gt(q') = 
G^iq') = 

G^iq') = 



Go 



P2 



+p, 



TO^,) 



(24) 
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Equating the single and double pole terms gives 
two relations. Using these and taking the double 
limit, and rUg 0, gives 



2MpGf{q^ = 0) = -ai + Nf-^F 



(25) 



2m„ 



9 1 9 



= r:2- -2- -Pi +^/:^(^- 
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The term proportional to Nfag/2'iT diverges in 
the chiral limit and there is no obvious way of ex- 
tracting the physical answer from it alone. Thus 
the method fails in the quenched theory. 

In the full theory, there are no double poles and 
an analogous analysis gives 



2MpGf{q^ = 0) 



Nf^iF- 



ml,- 



(26) 



which justifies the use of the anomaly relation. 

9. MASSES OF LIGHT QUARKS 

In order to extract light quark masses from lat- 
tice simulations we use an ansatz for the chiral 
behavior of hadron masses. Theoretically, the 
best defined procedure is xPT which relates the 
masses of pseudoscalar mesons to m^, m^, Ws- 

The overall scale /i in the mass term of Eq.l im- 
plies that only ratios of quark masses can be de- 
termined using xPT. The predictions from xPT 
for the two independent ratios are [6] [32] 

Lowest order Next order 
{mu + md)/2m, — — 

{md - mu)/ms 

In Lattice QCD it is traditional to make fits 
to the pseudoscalar spectrum assuming = 
j4^(mi + m2) and using either m.p or to set 
the scale. (The expression in Eq.5 is not relevant 
for this discussion since most quenched simula- 
tions have niq > »),.,/ 2.) A consequence of using 
just the linear term is that the ratio mg/m = 
25, i.e. these fits can be used to extract either 
m = (■/)),„ -|- m,i)/2 or /Ji., by using the physical 
masses for or mpc, but not both. (One would 
get a different number if 0{mg) and chiral log 
terms are included in the relation.) Furthermore, 
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since lattice calculations are done in the isospin 
limit, m„ = m^, therefore can be used to 

predict only one quark mass. The mass I prefer 
to extract, barring the complications of quenched 
xPT, is m as it avoids the question whether low- 
est order xPT is valid up to m,^. Akira Ukawa 
reviewed the status of m at LATTICE92 [28] and 
I present an update on it. 

To convert lattice results to the continuum MS 
scheme I use 

nicontiq*) = miau(a) [l--^{log(q*a)-Cm)](27) 

where the renormalization scale fi is the same 
as q* (defined in [29]) and chosen to be 7r/«, 
Cm = 2.159 for Wilson [30] and 6.536 for stag- 
gered fermions [31], and the rho mass is used to 
sot the scale. (I have not used the tadpole im- 
provement factor of Uq [29] in Cm and miatt as 
this factor cancels in perturbation theory and is 
a small effect otherwise.) The value of boosted 
I use in Eq.27 is [29] 



0.025 



(28) 



which is consistent with the continuum MS 
scheme value at Q = -n ja 



1 



g\Q) 

le^r^ /JoLn(f^) 



/3iLn[Ln<'^ 



A2 



(29) 



provided I use A = 245 MeV and 190Mey for 
n/ = and 2 theories respectively. Note that the 
choice of A, q*, and the constant 0.025 in Eq.28 
are interrelated and, at this order, one can trade 
changes between them. Finally, all the results are 
run down to Q = 2 GeV using 



m{Q) 
m[q*) 



5'(Q)- 5^(9*) /7i/3o- 70/3: 



167r2 



) .(30) 



The status of calculations of m(2 GeV) for 
quenched Wilson [35] [20] [36] [37] [12] [38], 
quenched staggered [39] [9] [40] [7] [8], dynam- 
ical irif = 2) Wilson [41] [42], and dynamical 
(nf = 2) staggered fermions [40] [43] [44] is shown 
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Figure 7. m extracted using data with the 
scale set by rrip. 
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in Fig. 7. I liave suppressed error bars as I want 
to first emphasize key qualitative features. The 
quenched staggered and the rif = 2 Wilson and 
staggered give m = 2 — 3 MeV and are roughly 
consistent; however, the quenched Wilson results 
seem to approach that value from above and even 
at /3 = 6.4 are significantly higher. (The recent 
result = 128(18) MeV by AUton et al. [45] for 
both Wilson and Sheikholeslami-Wohlert actions 
at /3 = 6.0 and 6.2 is consistent with the results in 
Fig. 7 once one notes that m = mg/25.) I believe 
that, at this stage, it is important to understand 
why the quenched results with the Wilson action 
are so different from the rest! 



An alternative to using 



Zg ^m^ to calculate the quark masses with Wil- 
son fermions is to use the Ward identity [34] [15] 

Za {A4t)P{0)) 



Zp 2 <P(t)P(0)) 



(31) 



Using the perturbative values for Za and Zp 
(with q* = Tr/a and boosted g'^ defined in 
Eq.28) the LANL Wilson data [12] gives m = 
3.53(10) MeV in contrast to m = 5.15(15) MeV 
shown in Fig. 7. The statistical errors are calcu- 
lated using a single elimination JK with a sample 



of 100 lattices of size 32^ x 64, so the difference 
is significant. The Rome collaboration [45] has 
found a similar discrepancy and argue that it can 
be resolved if one uses the non-perturbative value 
for Zp, which they advocate calculating using 
matrix elements of the operators between quark 
states in a fixed (Landau) gauge. Their results 
indicate that perturbation theory (including tad- 
pole improvement) fails for Zp. The two methods 
for extracting m give consistent results once the 
non-perturbative value of Zp is used. 

Having fixed m one can extract m^, mc, and mi, 
using, for example, K*, D, and B meson masses 
provided it is assumed that these masses are lin- 
ear in the light quark mass and in the heavier 
quark mass around the physical value. Alter- 
nately, one can use ■/Ji,/,. ,//(/'■ and T spectrum to 
get these quark masses directly without needing 
to extrapolate in the light quark mass. The re- 
sults for rric and rrib have been reviewed by Sloan 
[33] at this conference so I will only analyze the 
data for and compare these estimates to 25m. 
Note that the same data used to compile Fig. 7 
is used to calculate rrig from mx' and m^. The 
procedure for translating the value to 2 GeV in 
the MS scheme is also the same. The results in 
Fig. 8 show that that the estimate of from 
is systematically higher by 15 — 20% compared to 
25m. 

I will use the LANL data[12] to show that the 
systematic errors due to choice of hadron used 
to set the scale of the strange quark are now a 
dominant source of error. We find that, in the 
WS scheme at 2 GeV, = 25m = 129(4) MeV 
using Mfc, rrig = 151(15) MeV using Mfc, and 
■m, = 157(13) MeV using M^. Note that the 
latter two estimates give mgjni <^ 30, which is 
much closer to the "Next Order" prediction of 
XPT. The larger errors in these cases refiect the 
fact that on the lattice masses of pseudoscalar 
mesons are measured with much better statistical 
accuracy than those of vector mesons. 

10. CONCLUSIONS AND COMING AT- 
TRACTIONS 

The analysis of various quenched quantities 
show that the parameter 8 characterizing the 
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Figure 8. Comparison of extracted using 
and = 25m. The data are for quenched Wil- 
son simulations. 
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hairpin vertex in the rj' propagator lies in the 
range 0.1 — 0.2. As a result, for rrig < ma/2 
I find significant deviations from the lowest or- 
der chiral behavior in m^/m^. Therefore, I con- 
clude that extrapolation of quenched data, ob- 
tained with nig < ma/2, to the chiral limit can- 
not be done simply using full QCD formulae for 
quantities which have large contributions from 
enhanced logs. For quantities like the matrix cle- 
ment of the singlet axial vector current using the 
Adler-Bell-Jackiw anomaly, the quenched approx- 
imation fails altogether. 

The alarmists are busy calculating 1-loop cor- 
rections to other quantities to determine what can 
be extracted reliably from quenched simulations. 
Bernard and Golterman have extended the results 
presented at LATTICE93 [46] and calculated chi- 
ral corrections to the energy of two pions in a fi- 
nite box as derived by Liischer [47]. They find 
terms at 0(1) and 0(1/ L^), whose contribution 
could be substantial, in addition to modifications 
of the 0(1/ L^) term which is related to the tt — tt 
scattering amplitude [48]. Sharpe and Labrenz 
have extended the analysis of baryons to include 
the A decuplet [49]. Booth [50] and Zhang and 
Sharpe [18] have calculated corrections to heavy- 
light meson properties like and Bb- These new 



results and more data should provide a clearer 
picture of what is possible with quenched QCD 
by LATTICE 95. 

In the calculations of light quark masses we 
need to understand the factor of two difference 
between the quenched Wilson and staggered data. 
On the other hand, the quenched staggered data 
is consistent with the Uf = 2 Wilson and stag- 
gered data. The analysis presented here leaves 
open the question — is the agreement between 
quenched Wilson (and 0(a) improved SW action) 
data with the phenomenologically favored esti- 
mates of m (or equivalently m,,) fortuitous and 
an artifact of strong coupling? If so, then the 
Uf = 0. 2 staggered and nf = 2 Wilson data give 
an estimate of m that is 2 — 3 times smaller than 
the commonly accepted plienomenological value. 

The systematic errors due to the choice of 
hadron mass used in determining are signif- 
icant. Using mx* or m^ to extract m^ gives 
a <^ 20% larger value than that obtained from 
mjf . Even though the statistical errors are larger 
when extracting m^ from vector mesons, these 
estimates provide information beyond the lowest 
order xPT result m,, = 25m. Phenomenologi- 
cal estimates involving extrapolation to strange 
quark mass need to take this systematic differ- 
ence into account. 
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